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Abstract 

The application of the optimized expansion for the quantum-mechanical 
propagation in the anharmonic potential Ax 4 is discussed for real and 
imaginary time. The first order results in the imaginary time formal- 
ism provide approximations to the free energy and particle density 
which agree well with the exact results in the whole range of temper- 
atures. 
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1. The quantum mechanical Hamiltonian 



provides the simplest representation of anharmonic effects and is widely used 
in solid state physics, quantum chemistry and even in paraxial optics (where 
the Helmholtz equation reduces to the Schrodinger equation). Although the 
anharmonic oscillator (AO) cannot be solved exactly, the energy spectrum, 
as well as the time evolution amplitude can be calculated numerically to an 
arbitrary accuracy. For systems with many degrees of freedom such calcu- 
lation become unpractically time-consumming, especially in the case of the 
time evolution, so that approximation methods are desirable. 

A very promissing and simple approximation method is generated by the 
optimized expansion (OE). The method has been formulated for the effective 
potential in quantum field theory in the space-time of arbitrary dimension [|l|] . 
Since the AO is equivalent to the theory of a scalar field in the space-time of 
one-dimension (time) with a classical action given by 



A[x] 



T 



i(x 2 (t) -mV(t)) - \x\t) 



dt, (2) 



the OE for the effective potential can be used to generate a systematic ap- 
proximation scheme for the free energy of the system [[|, in the following 
called the optimized expansion for the free energy (OEF). It has been shown 
that the lowest orders of the OEF provide good approximations and the 
series converges to the exact free energy [EJ. 



2. Here we generalize the method to describe local properties of the system, 
applying the OE to the evolution amplitude. The real time propagator can 
be represented by 

rx(T)=x b 

(x b ,T\x a ,0) = / Dxe lA ^ x \ (3) 

Jx{0)=x a 

where the integral is taken over the functions which begin at x(0) = x a and 
end at x(T) = Xb- 

The most popular perturbation method is generated by expanding the 
propagator (0) into a series in the coupling constant A and performing the 
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path integrals analytically, order by order. Unfortunately the perturbation 
series for energy eigenvalues of the AO is asymptotic [§J, keeping thus a few 
lowest terms in the expansion provides reasonable approximation for very 
small values of A and short propagation times T only. Since the OE for the 
free energy shows much better convergence properties than the perturbation 
method, we expect also the approximations to the propagator generated in 
the OE to have much broader region of applicability. 

The OE consists in modifying the classical action to the form 



A 



mod 



X 



T r 



(x 2 (t) - uj 2 x 2 (t)) + e Q(c^ 2 - m 2 )x 2 (t) - Xx\t) 



dt, (4) 



where the unperturbed part corresponds to the harmonic oscillator with an 
arbitrary frequency, uj. The formal expansion parameter, e, has been intro- 
duced in such a way that for e = 1 the dependence on uj cancels and the 
modified action becomes equal to the classical one (Q); therefore, after cal- 
culating the quantity of interest to the given order in e we set e = 1. The 
exact result, obtained as a sum of an infinite series, does not depend on an 
arbitrary frequency, but a finite order truncation does. The advantage of 
the freedom can be taken to optimize the expansion: we make the nth-order 
approximant as insensitive as possible to small variation of u, choosing the 
value of uj which renders the approximant stationary. The optimal frequency 
changes from order to order, which improves convergence properties of the 
expansion scheme. 

The OE for the time evolution can be generated by expanding the am- 
plitude (H) into a series 



(x b ,T\x a ,0) 
r(x b ,T) 

J{Xa,0) 



(x b ,T) 

(Xafi) 
T ± 2 ( 



Dx e 



i ^[\{x 2 (t)~uj 2 x 2 {t))+e{\{u 2 -m 2 )x 2 {t)~\x i {t))]dt 



Dxe'Jo 



-a) 



1 + ie 



to 



iuj (x^-^-x 2 ) cos ujT — 2x a Xfo 




2 2 

uj — m 



1 -He- 



ar 



x 2 (t) 
2 

2 , T 



(2iri sinujT) 

-teX f T \£ 4 (t) + QiC 2 (t)1C(t) - ZK?(t)}dt + 0(e 2 
J o 







- Ax 4 (t)J dt + 0(e 2 ) 
\C 2 {t) +iJC{t)]dt 



(5) 
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where 

m = ^™t + x b sinuj(T-t) and m ^^^(T-t) 
sin ujT uj sin ujT 

Representing the evolution amplitude by 

(x b} T\x a} 0) = e m ^ T \ (7) 

we calculate W to n-th order in e and upon setting e = 1 we optimize the 
approximant, W n (x b , x a ,T), choosing uj to fulfill 

SW n (x a ,x b ,T) 

to " °- (8) 

The optimal choice of to in each order calculation, ensures that the propaga- 
tor about which we expand captures essential features of the system under 
investigation. 

The first order approximation is calculated to be given by 

.m 2 — uj 2 r T 



W^(x a ,x b ,T) = W°(x a ,x b ,T)-i / [C 2 {t)+iK,{t)]dt 

I Jo 

- iX f \C\t) + QiC 2 (t)K(t) - 3JC 2 {t)]dt, (9) 
Jo 



where 



, n , 1 / uj \ iuj \{x 2 a + x 2 ) cosuT - 2x a x b ] 



W u (x a , x b , T) = - In , , + lv a b > ^ (10) 

2 \2m sm col / Zsmul 

corresponds to a harmonic oscillator with a frequency uj. The optimization 
condition (f?) reduces to 



m 2 — UJ 2 rT 



[C 2 (t)+i}C(t)}dt+\ / T [£ 4 (t)+6i£ 2 (t)/C(t)-3/C 2 (t)]rft = (11) 
Jo 



2 

because of 

5W [> 1 r T 



( [£ 2 {t)+i!C{t)]dt. (12) 
Jo 



5u 2 2 Jo 

For the harmonic oscillator with a frequency m (A = 0) the optimization 
condition (|TT]) is solved by uj = m, so that the exact propagator is recovered 
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in the discussed approximation. Numerical results for the AO propagator 
(A 7^ 0) can be easily obtained, calculating however the integral for wave 
function evolution 



ip(x b ,T) = J (x b ,T\x a ,0)ip(x a ,0)dx a 



(13) 



the problem of highly oscillatory behavior is encountered. We decided there- 
fore to discuss first the quality of the approximation for the imaginary time 
propagator. 



3. The imaginary time evolution amplitude 



(x b ,P\x a ,0) 



jW(x b ,X a ,f3) 



/ Dxe~ A[x \ 



l(x a ,0) 



where t — it and the Euclidean action is given by 



A[x\ = / 
Jo 



-(x 2 (r)+m¥(r)) + Ai 4 (r) 
2 



dr, 



(14) 



(15) 



describes equilibrium properties of the system at temperature (5 1 . The trace 
of the imaginary time propagator defines the partition function 



/P fixa, /3) P 

dx a (x a ,(3\x a ,0)= JdxaJ ' Dxe~ A W = ]dx a e^\ 



(16) 



where W(x a ,f3) = W(x a , x a , (3) and the free energy can be obtained as 
= —YaZp/fl. The density matrix can be expressed as 



p(x a ,x b ) =Z p l / Dxe~ A \ x \ 



(17) 



and the average particle density is given by its diagonal element, p(x a 

P\X a , X a ) . 

The OE is generated by modifying the Euclidean action to the form 



j^mod 



X 



f 

Jo 



i(i 2 (r)+^V(r)) + e 



m z -u 2 )x 2 (t) + \x*(t) 



dr. (18) 



5 



The first order approximation for the imaginary time amplitude (0), ob- 
tained by an analytic continuation of the real time result (P), reads 



WW{x ai x h ,P) 



In 



2itsmhujf3 



+ 



uj [(x 2 + xl) coshcj/^ — 2x a Xb] 
2 sinhtu/? 



UJ 



[L 2 (t) + K(t)]dt - A / \L\t) + QL (t)K(t) + 3K{t) 2 ]dt (19) 



where 



r/ , x a sinhcut + Xb smh. uj(/3 — t) , . sinhtutsinhtuf/? — t) 

L ( t ) = and K{t) = - 



sinha;/3 



u sinho;/3 



(20) 

If we calculate the free energy by expanding the subintegral expression in 
Eq. [16| to the given order in e and performing the Gaussian integrals over x a 
we would obtain the series for the free energy 



m 2 — uj 2 


rl 


1 " 


3A 
+ — 


rl 


1 ■ 


2uj 


,2 H 


eP" - 1. 




,2 H 


- 1. 



which coincides with OEF 
requiring 



+0(e 2 
(21) 

The n-th order approximation is obtained by 



5£l 
Sou 



0. 



(22) 



which determines uj as a function of /3. The OEF provides information on 
global properties of the system only To discuss local properties we pro- 
ceed differently in this paper: we perform the optimized expansion for the 
propagator (OEP), imposing a local optimization condition 



6WW{x ai x b ,P) 



5uj 







(23) 



which determines uj as a function of /3,x a and Xb- This approach yields an 
approximations to the density matrix fllTD in a natural way. The given order 
approximation to the partition function is obtained by performing integration 
over x a in Eq. [16] numerically and the free energy is derived afterwards. 

Here we discuss approximations to the free energy and to the particle 
density for the one-dimensional AO, generated in the first order of the OE. 
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In Fig. 1 we show the results for the quartic oscillator (m = 0, A = 1): the 
free energy OEP (obtained by optimization of the imaginary time amplitude) 
is compared with OEF (obtained by optimization of the free energy) and the 
exact result calculated numerically. In the limit of high temperature both 
OEP and OEF approach the exact result, at zero temperature the approxi- 
mations coincide also, but the accuracy is the worst. At finite temperature 
OEP appears slightly better than OEF, differences between approximations 
are not large oving to the fact that only a small region of x a contributes 
to the partition function and it does not make a big difference whether the 
integration over x a is performed before or after optimization. The quartic 
oscillator is on the border between the single well (m 2 > 0) and double well 
(m 2 < 0) AO. Since the accuracy improves for increasing jjbs, the approxi- 
mations to the free energy for m 2 > are better than in the case of quartic 
oscillator. For the double well oscillator the accuracy is worse, but even in 
this case both OEF and OEP are satisfactory provided the wells are not 
very deep. The OEP has a great advantage of yielding directly the approx- 
imations to the density matrix ( |TTD and the particle density. The particle 
densities, calculated in the first order of the OE for the single well (m 2 = 1 
and A = 10) and double well potential (m 2 = —1 and A = .1) are shown in 
Fig. 2 and Fig.3, respectively. They are in good agreement with the exact 
densities calculated from the Schrodinger wave functions. 

The OE for the imaginary time propagation amplitude bears some simi- 
larities with the Feynman-Kleinert (FK) variational method [|5], ||, extended 
to a systematic variational perturbation theory for the free energy |7j . In the 
FK approach the exact partition function is expressed as 

dx -(3V cl (x ,l3) 



'(27T/3) 

with the classical effective potential, V c i(xq,P), defined by 

e -pv cl (x ,/3) = f Dx Jl 2 iTf3)5(x - x)e~ A[x \ (25) 



where xq = x = JjtZ^l^ j n the variational perturbation theory the same 



modified action fll8|) is used, but the partition function is obtained by per- 
forming numerical integration over xq in (p4[ ) with V c i(xq, j3) calculated to 
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the given order in e, i.e., the OE is applied to the classical effective poten- 
tial, while in our approach the partition function is expressed by W(x a ,(3) 
and integration over x a in ([II]) is performed. The free energy for the quar- 
tic oscillator obtained in the first order of the FK approach is compared in 
Fig.l with our results. In the limits of high and low temperature the FK 
results coincide with OEP, for intermediate temperatures the former agree 
better with the exact free energy than the later. To any finite order in e the 
results of the methods are different, the difference comes from the fact that 
V c i(xo, j3) is a function of the mean value of the coordinate, xo, and W(x a , (3) 
is a function of a starting point on periodic trajectory, x a = Xb = x(t = 0). 
The OEP can be thus directly applied to calculate the particle density which 
is a function of x a . In the FK approach the relation between xq and x d can 
be taken into account to calculate the particle density ||, but this requires 
an additionary numerical integration. The results of the FK method for par- 
ticle densities which are shown in Fig. 2 and 3 are of similar accuracy as that 
obtained in the OEP, the later are even better for the single well oscillator. 
This is remarkable since the FK approach requires two numerical integra- 
tions, while in the OEP only the integral for the partition function has to be 
performed numerically, which is an important advantage of our method in 
view of further applications to the systems with many degrees of freedom. 

It is worthwile to note that the FK method cannot be extended to the 
case of non-periodic trajectories in order to obtain non-diagonal terms of the 
density matrix and the propagation amplitude in the real time formalizm. 
The OEP ofers such a possibility in a natural way, we shall present a de- 
tailed discussion of the real time propagation in a separate paper. Both in 
the imaginary and in the real time formalism, the OEP gives a possibility 
to calculate corrections to the propagation amplitude, improving generated 
approximations in a systematic way. In this way the convergence properties 
of the approximation scheme can be estimated which is of great importance 
for the systems where the exact result is difficult to obtain. 
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Figure captions 



Figure 1. The free energy F of the quartic oscillator, (m 2 = 0, A = 1) in the 
first order of the OE, obtained by optimization of the imaginary time prop- 
agator (OEP, solid line), and of the free energy (OEF, dashed- dotted line) 
compared with the FK approximation (dashed line) and the EXACT result 
(dotted line), plotted vs. the inverse temperature (3. 

Figure 2. The particle distribution of the single well oscillator (m 2 = 1, A = 
10), obtained in the OE for the imaginary time propagator (OEP, solid line), 
compared with the FK approximation (dashed line) and the EXACT result 
(dotted line) at (5 — .1 and (5 — 5. 

Figure 3. Same as in Fig.2, but for the double-well oscillator (m 2 = — 1, A = .1) 

at = .25 and (3 = 5. 
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